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We study the gluon and ghost propagators of SU(2) lattice Landau gauge theory and find their 
low-momentum behavior being sensitive to the lowest non-trivial eigenvalue (Ai) of the Faddeev- 
Popov operator If the gauge-fixing favors Gribov copies with small (large) values for X\ both the 
ghost dressing function and the gluon propagator get enhanced (suppressed) at low momentum. 
For larger momenta no dependence on Gribov copies is seen. We compare our lattice data to the 
corresponding (decoupling) solutions from the DSE/FRGE study of Fischer, Maas and Pawlowski 
[Annals Phys. 324 (2009) 2408] and find qualitatively good agreement. 
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1. Introduction 

Studies of the elementary two and three-point functions of Landau-gauge Yang-Mills theory 
have always been an interesting topic on the lattice. In particular the gluon and ghost propagators 
became one in 2005 when it was apparent that (continuum) functional methods [1] and lattice ap- 
proaches [2] disagree in their findings for the propagator's low-momentum dependence. Since then 
many efforts have been made, both on the lattice and in the continuum, to verify and understand 
this discrepancy [3-5]. 

A possible solution was proposed by Fischer, Maas and Pawlowski in 2008 [4]. They investi- 
gated the (truncated) Dyson-Schwinger equations (DSEs) of the gluon and ghost propagators, and 
in addition also the corresponding functional renormalization group equations (FRGEs), and found 
there is not a unique solution to the system of equations but a one-parameter family of decoupling 
solutions. A particular solution is chosen by the value set for the ghost dressing function at zero 
momentum, /(O). In the limit /(O) — )■ °°, one obtains the scaling solution that was found before 
[1], but not yet on the lattice. 

Although this proposal may be attractive to understand the discrepancy, lattice studies have not 
delivered evidence. In this contribution we show that a part of this family of (decoupling) solutions 
can be reproduced on the lattice, at least qualitatively, as far as possible on a finite and rather coarse 
lattice, and to the extent computational resources allow. We are also limited by our approach which 
allows only for mild variations of the ghost dressing function at low momenta. For this variation we 
utilize the lowest non-trivial eigenvalue Ai > of the Faddeev-Popov (FP) operator and show that 
the low-momentum behavior of both the ghost dressing function and the gluon propagator changes 
with the average Ai of the selected gauge-fixed (Gribov) copies.^ Interestingly, these changes look 
qualitatively the same as for the corresponding subset of decoupling solutions. 

2. Setup 

We restrict our study to SU(2) lattice gauge theory (Wilson plaquette action), and also consider 
only one lattice size (56^*) and gauge coupling (j3 = 2.3). This is fully sufficient for our purposes 
and allows us to analyze (with reasonable computing time) data for the gluon propagator and the 
ghost dressing function where their momentum dependence starts to become flat, and this for a 
large enough number of Gribov copies such that a correlation between Ai and the ghost and gluon 
propagators can be seen. Moreover, by restricting to one lattice spacing and volume no further 
effects (finite volume, renormalization, discretization) interfere. 

Our results are for an ensemble of 80 thermalized gauge field configurations. These are sepa- 
rated by 2000 thermalization steps, each involving four over-relaxation and one heatbath step. This 
number turns out to be sufficient as no apparent autocorrelation effects are seen in the data^. For 
every gauge configuration there are at least A'copy = 210 gauge-fixed (Gribov) copies, generated 



'Alternatively, one could also directly constrain J(p) at some p > to select Gribov copies (see [5]). 
At the conference results were presented for 60 configurations. To improve data and also to verify that our results 
do not suffer from large autocorrelation effects, another (independent) chain of 20 configurations was generated. The 
results on either chain are fully compatible with each other. A binning analysis is applied to estimate statistical errors 
nonetheless. 
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using an optimally-tuned over-relaxation algorithm. To ensure that these copies are all distinct, the 
gauge-fixing always started from a random point on the gauge orbit. Only a few Gribov copies 
were found twice. 

For each Gribov copy we calculate the lowest three (non-trivial) eigenvalues < Ai < A2 < A3 
of the FP operator and use X\ to classify copies. For each gauge configuration, the Gribov copy 
with the lowest value for X\ is labeled lowest copy (£c), while that with the highest X\ is called 
highest copy {he). The first generated copy, irrespective of Ai, gets the label ^r^f copy (fc). It 
represents an arbitrary (random) Gribov copy of a configuration. To compare with former lattice 
studies on the problem of Gribov copies we also consider best copies, i.e., those copies with the 
best (largest) gauge functional value 

Fu[8] = ^'L'L^^Trg,U.-,gl^ (2.1) 

^'^ X 11=1 

for a particular gauge configuration U = {Uxfi}. Here g = {gx} denotes one of the many gauge 
transformation fields fixing U to Landau gauge. 

The gluon (D) and ghost propagators (G) are analyzed separately on those four sets of Gribov 
copies. We apply standard recipes for their calculation. For the ghost propagator, though, we will 
primarily analyze its dressing function / = p^G. It parametrizes the deviations from the tree-level 
(infrared diverging) propagator and is thus better suited for our purposes. When quoting momenta 
in physical units we adopt the usual definition a/7^(A;^) =2sin (7rA;^/L^) [withA;^ G (— L^/2,L^/2] 
and L^ = 56], assume for the string tension ^/a = 440 MeV and use aa^ = 0.145 for j8 = 2.3 from 
Ref. [6], where a denotes the lattice spacing. 

3. Results 

Data for the ghost dressing function and the gluon propagator is shown in Fig. 1. Looking 
there first at the left panel, one clearly sees the choice of Gribov copies affects the momentum 
dependence of the ghost dressing function at momenta p^ < 0.2 GeV^. Points for the different sets 
of Gribov copies deviate systematically and the effect increases the lower the momentum. For the 
ic data we find the strongest enhancement for the ghost dressing function towards zero momentum. 
The he data shows the weakest enhancement, and this data also almost coincides with the be data 
(see also the discussion below). On the other hand, for momenta above 1 GeV no Gribov-copy 
effects are seen. 

Interestingly, also for the gluon propagator we see Gribov-copy effects at the lowest momenta 
(right panel of Fig. 1). Due to the larger statistical uncertainties (as typical for this propagator), 
these effects are less significant however. A much enhanced statistics would be desirable, but 
unfortunately this is beyond our cun^ent resources of computing time. Nonetheless, a systematic 
deviation of the £e (red) and he data (green) is seen for momenta p^ < 0.1 GeV^, while for larger 
momenta the data points for all sets agree within errors. We also see that the be data is suppressed 
compared to the fe data, in agreement with earlier studies [3, 7]. 

The coiTclation between the propagator's low-momentum behavior and Ai is even better seen 
when looking at the values for Ai and the corresponding ("measured") values for the propagators, 
that is sepai^ately for each Gribov copy. Such scatter plots are shown in Fig. 2; the top panel shows 
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Figure 1: Ghost dressing function (left) and gluon propagator (right) versus a^p^. Full squares (open circles, 
crosses, full triangles) refer to tc (fc, be, he) data; shown is the raw lattice data, that is no renormalization 
has been applied. A zoomed-in plot for the gluon propagator improves the visibility of the low-momentum 
region. For the same reason, points at p = are slightly shifted. Physical momenta are given at the top. 



the ghost dressing function at the lowest finite momentum {a'^p^ ~ 0.01258) and the other two the 
gluon propagator at a-^p'l (middle) and at p^ = (bottom). To increase visibility of the Ai depen- 
dence, we also show averaged values (colored bars) over adequately chosen, partly overlapping Ai 
intervals. The bar length equals the Ai interval and the bar height reflects the statistical uncertainty 
of each average (marked by a line). 

For the ghost dressing function we clearly see the data points to grow if Ai is decreased and 
vice versa. In particular towards small Ai the effect becomes large. We also find (middle panel) 
the fc points to lie above the be points, as expected from other studies, but this dependence on the 
gauge-functional value adds to the dependence on Ai . This may explain why the be and he ghost 
dressing functions (accidentally) coincide as noticed above. 

In comparison, the gluon propagator data comes with much larger statistical fluctuations which 
makes it hard to draw finite conclusions. Nonetheless, a trend is seen in the data: for large Ai the 
data points tend to lower values than for smaller Ai . Though, it is hard to decide if this trend persists 
for very small Ai . From the middle and bottom panels of Fig. 2 one could conclude either. Note 
that this ambivalence is also reflected in the gluon propagator data shown in Fig. 1. There we see 
the order of the ie and fe points at p = and low p > being partly inverted. Though, given the 
findings below (Fig. 3), this may be just a statistical artifact. 

Let us remind that our averages for D{0) result from about 80 individual "measurements", 
one per gauge copy, while each Monte-Carlo history point for D{p > 0) itself is an average over 
all possible momentum directions with same a^p^. For the lowest finite p^ there are already 4 
directions one averages over. This explains the larger statistical noise for D{0). 

To cross-check if the Ai -dependence we see for the gluon propagator is not just a statistical 
artifact, we perform additional calculations of the gluon propagator on all Gribov copies with the 
second lowest and second highest Ai . Those additional sets of (2 x 80) Gribov copies (labeled s£e 
and she in what follows) are distinct from the sets of lowest and highest copies (ie and he) analyzed 
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Figure 2: Ghost dressing function (top) and gluon propagator at the lowest finite momentum p\ (middle) 
and zero momentum (bottom) versus Aj . In the background we show the "measured" values, separately for 
each Gribov copy (scatter plot), and in the foreground averages (colored bars). The middle line of each bar 
marks the average over the shown Ai -interval (bar length); the height reflects the statistical uncertainty. Note 
the different scales of the ordinates, and the partly overlapping ranges of Aj for the three panels showing 
data for our different sets Gribov copies (first, lowest, best and highest copies: fc, £c, be and he). 
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above, and if there is a dependence on X\ , one should also see it when comparing sic and she 
data. And in fact, also this data clearly exhibits a Ai -dependence at low momenta (see Fig. 3). The 
combined ic and sic data and the combined he and she data (colored en^or bands in Fig. 3) currently 
give the best impression of this dependence. Note that such a combination of data is justified, as 
there are no coiTclations visible between data from different copies of the same configuration, 
and by construction these sets of Gribov copies come also with similar small or large values for 
Ai: The averaged Ai values on these four sets of Gribov copies are: (Ai)^ = 1.43(9) x 10^^, 
(Ai),fc = 2.14(9) X 10-"^, (Ai)/,c = 20.33(6) x lO""* and (Ai),/,^ = 19.89(4) x 10""^ (lattice units). 
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4. Conclusion 

Our lattice study shows that the low- 
momentum behavior of the Landau-gauge 
gluon and ghost propagators can be changed 
(slightly, though systematically) on the lat- 
tice by constraining the lowest non-trivial 
eigenvalue Ai of the FP operator. If we re- 
strict Ai to be small (large) for each Gri- 
bov copy, the ghost propagator at low mo- 
mentum gets enhanced (suppressed), while 
it is not at all affected at intermediate or 
large momentum. Interestingly, also the 
gluon propagator (D) is affected at low 
momentum, but in comparison to the ghost 
dressing function the effect is smaller. Cur- 
rently the effect is best seen if one com- 
bines the data from Gribov copies with 
lowest and second lowest Ai, and from 
copies with highest and second highest Ai 
as shown, e.g., in Figs. 3 or 4. 

In Fig. 4 we also confront our data for the ghost dressing function and the gluon propagator 
with a coiTcsponding pair of decoupling solutions from [4]. These DSE solutions are approximately 
those where the boundary condition on the ghost dressing function was set to /(O) = 3.4 and 
/(O) =3.8, respectively. For the comparison our data has been renormalized relatively to the given 
decoupling solution, by applying a common renormalization factor (Zj and Z^) to the respective 
data. Since truncation effects are expected to become important, these two factors were chosen 
such that the he-data, points (green triangles) agree with the /(O) = 3.4 curves (green) at the second 
lowest finite momenta (this point is highlighted by green circle in the figure). One could of course 
chose any other point, which would result in a similar figure. But at the moment our comparison 
is only qualitative anyway. Nonetheless, the surprisingly good agreement between the so different 
approaches to the gluon and ghost propagators of Landau-gauge Yang-Mills theory is encouraging. 

More details will became available soon in a revised version of [8]. 



Figure 3: Gluon propagator versus a^p^. Data is shown 
separately for the sets of Gribov copies with lowest (ic), 
second lowest (sic), highest (he) and second highest (she) 
value for Ai. Error bands in the background are for the 
combined ic and sic data (red) and the combined he and 
she data (green). At p = points are slightly shifted. 
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Figure 4: Ghost dressing function (top) and gluon propagator (bottom) versus p^. Full symbols refer to our 
lattice data and lines to four selected decoupling (DSE) solutions from [4]. Note that the order of the gluon 
propagator lines at low momenta changes somewhere between 7(0) = 3.8 and 10. 
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